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Big Whorls Carry little whorls 

By J.C.R. Hunt 1 , J.C. Buell 2 , and A. A. Wray 2 


The space-time structure of homogeneous isotropic turbulence has been studied us- 
ing a direct spectral simulation on a 64 3 mesh at a microscale Reynolds number 
of Re\ = 48. Steady body forces were applied to a few low- wavenumber modes 
to make the flow statistically stationary. The results for the two-point space-time 
correlations of velocity and pressure (R n and R pp ) show that the auto- correlations 
of u and p are positive and have about the same integral time scales, and that the 
spatial correlation R pp {r) ~ i?n(0, r, 0). The form of R pp (r ) and the result that 
p 2 /(pu q) 2 ~ 1.0 agree fairly well with the Hinze/Batchelor (Hinze 1975) results. 
The three-dimensional energy spectrum for small space-time scales of velocity and 
pressure are consistent with the concept of large eddies advecting the small scales 
with a random Gaussian velocity (rms value of one component is u 0 ): the wavenum- 
ber ( k ), frequency (u?) and the energy spectrum £(k,u>) are related to the spatial 
energy spectrum E(k) by the formula 

£(fc,u>) = E(k)exp[-\u 2 /(aku 0 ) 2 ]/(y/2ir(aku 0 )), 

where a ~ 0.4-0.5. The same form is found for the pressure spectrum. Extrapolating 
this result to high Reynolds number implies that in the inertial range the frequency 
spectra are <£n(u>) oc (eu 0 ) 2 / 3 u>~ 5 / 3 (Tennekes 1975), and <£ pp (u>) oc (eu 0 ) 4 ^ 3 u;~ 7 / 3 . 


1. Introduction 


1.1 Background to the project 

Recent reports on the structure of turbulence and its dynamics have tended to 
emphasize the representation and the dynamics of the spatial structure rather than 
the temporal structure of the turbulent velocity field. However, closer examination 
of dynamical theories (e.g., Kraichnan’s, described by Leslie 1981) show that they 
are always based on certain assumptions about the temporal structure, which have 
not been subjected to detailed examination or computation. 

The temporal structure of turbulence also needs to be understood in order to 
develop models for the effects of turbulence and fluctuating pressure fields on the 
generation of surface and internal waves (Carruthers & Hunt 1986) and on dispersion 
of pollutants, particles and bubbles. One of the ways of developing practical models 
for these problems is first to develop a stochastic representation of the velocity field. 
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FIGURE 1. Flow field composed of large scale eddies and small scale eddies that 
are transported by the large scales. The small eddies also deform on a larger time 
scale. 


The current research on the spatial structures (such as that of Moin 1987 and Moin, 
Adrian & Kim 1987) needs to be supplemented by information on how the turbulent 
velocity field evolves in time. (Some initial suggestions were put forward by Turfus 
&: Hunt 1986, but they will be superceded by the present work.) 

The clearest recent account of the temporal structure of turbulence is given by 
Tennekes (1975) and Tennekes & Lumley (1979). Essentially, the current under- 
standing is that the largest scales of turbulence with velocity and length scales uo 

and L are unsteady and change on a time scale L/u 0 . Meanwhile, smaller eddies on a 

1 / 2 

scale t with velocity u(C) ( defined by the structure function [|u(x) — u(x + r)| 2 ] , 

where |r| = t ), have two time scales: the Eulerian time scale 

te(() ~ e/u 0 , (1.1a) 

for the eddy to be advected by the large eddies past an observing point (moving 
with the mean flow), and the Lagrangian time scale 

r L (t) ~ t/u(£), (1.16) 

for the velocity field on a scale (£) to change as it is advected by the large scale 
eddies. See figure 1. 

Thus at a point (moving with the mean flow), the dominant time-dependent 
phenomenon (which determines du/dt) is the rapid, random movement of small 
scale eddies past the observer by the large eddies. On the other hand, for a point 
moving with the fluid, the only change of velocity is the slow change caused by 
the dynamical processes at the scale l. These differences are best, defined by the 
Eulerian microscale 

- 1/2 


T<®> = \{du/dt) 2 / U 2 ] 


9 


(1.2a) 
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and the Lagrangian microscale, following a material element, 




= [(du/dt) 2 /u 2 ' 


^tT 1 / 2 


(1.26) 


For high Reynolds number turbulence these can be related to the basic time scale 
L/u 0 by 

r (B) ~ Re~ 1/2 L/u 0 , (1.3a) 

and 

r (i) ~ Re-'^L/uo, (1.36) 

where Re — u 0 L/v. Tennekes (1975) found that these estimates were consistent 
with the grid turbulence measurements of Shlien & Corrsin (1974). Another impor- 
tant consequence of Tennekes ’s analysis is that for high Reynolds number turbulence 
in the inertial subrange, the frequency spectrum <f> n(u>) of one component of the 
velocity, say «j, measured at a point (moving with the mean flow) has a form which 
is different from the universal Lagrangian form 


<t>ii\v)= C’Zeu) 2 , (1.4a) 

namely 

4>n = C'5(ei/ 0 ) 2/3 w _5/3 , (1.46) 

where is a constant for a given flow. Since this form depends on the large scales 
which do not have (even approximately) a universal character, it is not likely that 
cj is a universal constant, unlike the constant for the one-dimensional wavenumber 
spectrum <An(^i) or energy spectrum E(k). The prediction (1.4b) has not to our 
knowledge been properly tested experimentally or computationally. 

No detailed analyses have been developed for the four-dimensional wavenumber- 
frequency spectrum tensor ^ij(k,u;) or its equivalent energy spectrum £(k,u>), de- 
fined by 


*y(k,w) 


and 


1 

(2tt) 4 


III /^^(x + M + T )e— 

£(*, w)= J ¥«(k,u;)<L4(k). 


’ dr dr , 


|k|=* 


(1.5a) 

(1.56) 


Note that 



£(k,u))duj. 


(1.5c) 


However, Carruthers & Hunt (1986) made a rash speculation based on Tennekes’s 
time scale argument (in order to estimate internal wave generation) that 


- ^ij(k)6(w ± ku 0 ), 


where $,j(k) = f ^ij(k,w)dw. (1.6) 

J — OO 



80 


J.C.R. Hunt, J.C. Buell , and A. A. Wray 


Measurements have been made of one-dimensional space-time correlations R(r , r) 
of velocity and pressure fluctuations in grid turbulence and in boundary layers 
(Favre et al. 1956; Wills 1971), and in a pipe (Sabot & Comte-Bellot 1976). The 
general features have been described by formulae of the form (in coordinates moving 
with the mean flow) 

Rn{ri,r) ~ #n(r)F(r, r, u 0 , L) 

where F(r =0) = 1, and F decreases with r more slowly as r increases. See 
Hinze (1975 p. 416). Similar forms are observed for fluctuating pressure. Wills 
(1971) also took the Fourier transform of these measurements and calculated the 
one-dimensional wavenumber-frequency spectra 4> pp (ki , u;). When converted into a 
frame moving at the same mean speed as the travelling pressure field, he found that 
his result could be expressed as 


(i-7) 

where / ~ exp ( — \ tt ~2 ) , and uq is the rms value of the streamwise velocity com- 
ponent (at y/S ~ 0.5). Moser & Moin (1984) computed space-time correlations for 
channel flow. 

An interesting quantity that can be derived from these correlations is the Eu- 
lerian time scale (in a frame moving with the mean flow) of the i th velocity 

component. [Since the measurements are not extensive enough for T (£) to be com- 
puted from integrals of Ru{t), T^ e ^ has to be estimated from the value of r at 
which i?u — 1/e.] The results of Favre and Moser & Moin indicate that the value 
of T\E\ normalized in terms of the rms velocity and the integral length scale L\ l in 
the flow direction has a range of values given by 

P m = Tl m ^/L", ( 1 . 8 ) 

where 1 < < 2. Snyder & Lumley (1971) inferred a value of — 3 from 

their measurements of small particles in grid turbulence. 

Using the same normalization for the Lagrangian integral time scale T^\ defined 
by 

P {L) = Tl L) \fii/L[\ ( 1 - 9 ) 

experiments and simulations indicate that (3^ ~ 1 for grid turbulence (Snyder & 
Lumley 1971), and boundary layer turbulence (Durbin & Hunt 1980). A random 
Fourier mode representation for homogeneous turbulence has produced a value of 
P (L) ~ 1/3-1/2. 

So the current position is that the magnitudes of /3^ E ^ and are not well 
defined. 



Big Whorls Carry little whorls 


81 


1.2 Aims of the project 

The aim of the research presented here is to explore the space-time structure of 
homogeneous turbulence by computing and then interpreting the two-point spectra 
and correlations of the velocity and pressure fields. Many of these statistics are of 
considerable practical importance, as indicated in the previous section. In partic- 
ular it is of interest to compare the different time and length integral scales and 
microscales for Eulerian and Lagrangian quantities, and to compare the space and 
time spectra. 


2. Some Theoretical Ideas for Guiding Interpretation 

2.1 Velocity spectra 


2.1.1 Simple non-interacting eddies 

Consider a set of smooth eddies on a line, each located at a position x n , with 
a spatial scale £ n , a turnover time w" 1 , and a random phase e n . We consider, for 
simplicity, one component of turbulence, uj ; then a typical form for uj is 


«.(«,*)= f>n exp 

71—0 


(* ~ =Cn ) 

2P n 


21 


cos(u>„f + e n ), 


( 2 . 1 ) 


where a„ are random amplitude coefficients which are independent, a n a m = 6 nrn a £ . 
We assume the mean of ta n is and its random element has variance . 

The cross correlation at points x and x + r and times t and t + r, when averaged 
over a length X ~ ^2 £ n and over a sufficiently long time, is 

/ — AT 2 

#n(r,T)= 2Y5^«^nexp(-^-)exp(-i<r2 n r 2 )cos(u^r). (2.2a) 

n=0 n 

So in general the space-time correlation is not simply related to the spatial correla- 
tion Rn(r). But, for the case of just one eddy scale L , where N = 0, i2jj reduces 
to 

/— 2 

#n(»%T)= ^ a o Lex v(-{tf) exp(-i<T^ 0 T 2 )cos(u;oT), (2.26) 

or 

Rn(r,r) = /2 n (r,r = 0)/(t). 

Then the time and space correlations are separable. Similarly the wavenumber and 
frequency spectra are not linked, for this type of flow field. 

2.1.2 Small eddies moved by a large-scale velocity field 

Now consider a set of small eddies moved by a random large scale velocity U. 
These might be the small scale components of turbulent flow. Again, take one 
dimension: Let 

M N 

E E ft-nm exp [t(fc n *-|-u> m f)], 

— Af n=-~N 


U = 


(2.3) 
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where u; m = Uk m + u>' m , a n is again a set of random amplitude coefficients, and w' m 
is the frequency in the adverted frame. 

Let the probability density function of U be p{U). Then 

/ M N 

X X l a nml ex P [i(k n r + k n rU + u' m T)\p(U)dU. (2.4) 

m=-M n=-N 

It is interesting to consider two possibilities for p{U). 

(a) p{U) is Gaussian with variance U^: 


M N 

Rni r,T)= -£ E 


(2.5) 


m= — A4 n~ — N 

Then the one-dimensional wavenumber frequency spectrum is 
- Af JV -oo too 

«*>'-) = <555 E E / / e<-^ 

l 2 ** m= — M n= — N " J ~°° 


= £ ( “ 2 t ;“,| ^]/ (v^|t.H>o)- 


(2.6) 

For high wavenumbers, the advective time scale (k i Uo)~ i is much less than the 
time scale for the eddies to change as they move. Consequently |u>jj <C l^iCol- By 
taking large enough values of N to give a dense distribution of Fourier modes, a 
spatial spectrum <f>{ki) from ^(fci,u>) can be defined. From (2.6), 

V>(fci,u>) = <£(fci)exp^— / (V^lfcil^o). ( 2 - 7 ) 

Note that _ f00 

u 2 = <f>(ki)dki = / / xl>(ki,u;)dkidu. (2.8) 

J — oo J — ooJ — oo 

In a three-dimensional isotropic velocity field, u> n = U • k n + u)' n , and for a 
Gaussian distribution of the 3-D large-scale velocity U, it can be shown that the 
relation between the k-u energy spectrum and the k energy spectrum is 


£(k, u>) — E(k)ex p {-wui)h^ ku ^ 


(2.9) 


where Uq = ||U| 2 . 

(b) The large-scale velocity has a constant magnitude (forwards or backwards): 
In this case p(U) = \ [6{U - U 0 ) 4- S(U + U 0 )} , and this leads to 

fl„(r,r) = Y, \al\e ik - T e^-- r cos(rkM, 
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FIGURE 2. Similarity of the length scales of velocity and pressure fields induced 
by large vortices. 


and 

[^(u> — ktU 0 ) + + kjUo)]. (2.10) 

2.2 Pressure 

For large eddies having a form similar to vortices (such as analyzed by Townsend 
1976), the large scale distribution of pressure across the eddy, say in the X] direction, 
has a distribution similar to that of tt 2 (xj) (figure 2). Consequently, the pressure 
correlation R pp (r i) is expected to have the same scale and form as ^ 22(^1 )• 

The pressure fluctuations at small scales have been found in previous investi- 
gations (e.g. George et al. 1982) to be caused by the motions of eddies on that 
scale, and not (for example) by interactions between small scales and large scales. 
The former (which is observed) gives <j> pp {k j) a e 4 / 3 k 1 7 ^ 3 , while the latter gives 
ippiki) oc (tf 2 /I 2 )e 2 / 3 fc~ 11/3 . Thus it is expected that, as the large scale ed- 
dies advect the small eddies, they also advect the small scale pressure fluctuations. 
Therefore the k- u> power spectrum of pressure fluctuations is expected to have the 
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£ pp (k,u>) = E pp {k)ex p(- / (V2nkU 0 ). (2.11) 

This is close to the values measured by Wills (1971). 

3. Method 

To obtain the frequency spectra and time correlation of a turbulent velocity field, 
a long time record is necessary; a decaying turbulent field is not appropriate because 
the amplitude decay does not allow a long enough time record. It is necessary to 
generate a statistically stationary velocity field. 

The method used in these computations is to introduce a steady nonuniform force 
field F(x) at the largest scale of the flow (defined to be 2n in all three directions), 
which induces a mean flow at this scale. The mean flow is unstable, allowing 
instabilities to grow. This leads to a chaotic structure with motions at all scales 
down to the Kolmogorov microscale. If the initial conditions are chaotic, the force 
field will maintain the turbulence, but not necessarily at the original amplitude. 
Given a spatially periodic flow (whose maximum scale is fixed) any steady force 
distribution will eventually give rise to a stationary turbulent flow. This flow is 
determined by F and the viscosity u. 

The computations were performed on a N x x N y x N z mesh with periodic boundary 
conditions. The force spectrum j E f (k) is chosen to have contributions from values 
of k = y/2, i.e., k = (±1,0, ±1), (±1, ±1,0), and (0,±1,±1). To generate isotropic 
turbulence, all the moments of the body force F(x) would have to satisfy isotropy 
conditions. This could only be done approximately, by ensuring that moments up to 
second order of the Fourier coefficients of F, f, satisfied isotropy. We also specified 
that / 1 / 2/3 = 0. F was chosen to be solenoidal to avoid the generation of large 
pressure fluctuations. For each k such that |k| = y/2 the amplitude of the forcing 

was set so that yj /,(k)/*(k) = 0.2. The viscosity v was .02. 

The 3-D Navier-Stokes equations were solved under the above conditions using 
the spectral code of Rogallo (1981) modified to include the body force F. The 
primary results of the computation were the time-dependent, 3-D spatial Fourier 
coefficients of the velocity, u, and pressure, p, fields, defined by 

N*/2 N,/ 2 NJ 2 

Uj(x,<)= £ £ ^ «j(M)exp(ik-x), (3- 1 ) 

ki = -N x / 2 fcj = -N,/2 k s ~-Nj2 

where N x = N y = N z = 64. A similar definition applies to p(k,t). 

After storing «j(k,<) at M time intervals A t s , the individual time series are 
filtered with a cosine “window” within the period T — MAt s . Fourier transforms 
in time yield the 3-D space-time Fourier coefficients «j(k,u> m ), defined by 

M/2 

Uj(k,<)= ^2 «j(k,u; m )exp(iu) m <), 

m= — M/2 


(3.2) 
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where u> m = 2mrn/T. 

The time interval A t, used for the computation of frequency spectra and time 
correlations must be small enough to resolve the shortest time variation at a point. 
This minimum Eulerian time scale i s defined by the fast large eddies carrying, 
at a velocity u 0 , the smallest scale eddies £ro 1 past the point, i.e., ~ ^Koi/«o- 

(It is Re A ' times the magnitude of the Eulerian microscale!) In our computations 

At (based on CFL limitations) turned out to be less than O.lr^^, thus we took 
At, = 10At. 

The flow was initiated by suddenly imposing the force field on a given initial 
random isotropic velocity field u(x, t = 0). It was found that uj initially increased 
or decreased, but eventually oscillated with small amplitude about a stationary 
value; then the time series were collected to obtain the statistics. 

From the Fourier coefficients, space-time spectra and correlations may be com- 
puted, using the normalized length (27r) of the box. 

(i) The wavenumber-frequency (k,u>) energy spectrum tensor, 


*ij(k,u>) = tti(k,o>)u*(k,w), 

(3.3) 

and similarly for '® , pp (k,u>). 

(ii) The energy spectrum tensor, 


$ij(k) =Ui{k,t)u*(k,t). 

(3.4) 

(iii) The energy spectra (summing over spherical shells in wavenumber space; the 
equivalent integral expressions for a continuum of wavenumbers are given in (1.5)), 

fc- j<|k|<l+i 

(3.5a) 

E(k) = £ *,,(k). 

i-l<|k|<*+l 

(3.56) 

(iv) The one-dimensional space-time correlation, 


Rij (r,r) = ££*y< k -“> £ ‘ <k ' r+ "" r) - 

k m 

(3.6) 

(v) One-dimensional spectra, 


*2 *3 

(3.7) 


4. Results of the computations 

Table 1 gives most of the major statistics of the computed turbulent flow. Note 
that the Reynolds number based on the Taylor microscale, A, is only 48, so that we 



86 


J.C.R. Hunt, J.C. Buell, and A. A. Wray 


do not expect to see many features of high Reynolds number turbulence. However 
the Reynolds number is large enough so that the effects of advection of small scales 
can be investigated. Note also that the turbulent velocity field is not quite isotropic; 
the lack of isotropy in the moments ( u\/u\ ~ 1.1) for a given realization is of the 
same order as that usually found in wind tunnel experiments. This probably means 
that the forcing is not quite isotropic. The ratio of the integral scales L' u , L 22 , 
computed from the one-dimensional correlations -Rn(ri,0,0), ^ 22 (^i ? 0, 0), is 1.64. 
This is less than the value of 2.0 required by isotropy. By comparison the decaying 
strained turbulence computed on a 128 3 mesh by Rogallo (1981) and others was 
effectively defined to be isotropic at its initiation. 

Figure 3 shows the unfiltered and filtered time variations of typical large scale 
and small scale Fourier coefficients. Note the rather sudden changes in the small 
scale filtered velocity. The time scale t*. for such a change is consistent with large 
scales advecting the small scales (r*. ~ 2i r/(uok)). 

Figure 4a shows the contours of the energy spectrum £(k,us), and figure 4b the 
various statistics of the variation of £(k, u>) with u> for different values of k. The 
mean and skewness are close to zero. The standard deviation defined by 


as 


[/: 


u> 2 £{k,u>)du> 


//: 


1/2 


£(k,u ;)dw 


(4.1) 


is described by the formula &£ = aku 0 , where i/. 0 = y \u? . The flatness (or kurtosis) 
is close to 3. Consequently £(k, u>) can be approximated by 


£(fc,u>) = E(k)ex p 


2(aku 0 ) 2 


(\/ 27 r(afcu 0 )), 


(4.2) 


where a ~ 0.51 for kL\ < 10, and a ~ 0.4 for fcTj > 15 (figure 4c). A similar 
result for the k-u> energy spectrum for pressure has been found with a p ~ 0.51 for 
kLi < 10, and a p ~ 0.45 for kL\ > 15. This result for £ pp (fc,u>) is in approximate 
agreement with the space-time correlation measurements of Wills (1971 ) for pressure 
fluctuations on the wall below a boundary layer. 

An important consequence of (4.2) is that by assuming its validity at high 
Reynolds number we can calculate the Eulerian frequency spectra for the veloc- 
ity (f> u(u7) and the pressure <f> pp ( u>). Since <f> n(w) = f / £(k,u)dk and E(k) = 
afcc 2 / 3 fc -5 / 3 , where a* is the Kolmogorov constant, we have 


4 >\ i(^) 


2 2 / 3^5 


afc(e1io) 2 ^ 3u, 5 ^ 3 - 


Taking a = 0.51 yields 


^„(u;) = 0.085ajfe(cu 0 ) 2/3u;-5/r3 - 


(4.3) 
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FIGURE 3. Typical time variation of computed Fourier coefficients over the period 
of time used for frequency calculations. The real part of Ui(k x ,k y ,k z ) is shown. 
The data is filtered to allow processing by discrete Fourier transform methods: the 
filtering is most apparent at the interval boundaries, (a) data at low wavenumber, 

tii (1,1,1): unfiltered, filtered; (b) filtered and unfiltered data at high 

wavenumber, tt2(4,4, 19). 


Similarly if E pp (k) = afc p e 2 / 3 A: 7 / 3 (George et al. 1984), we deduce that 
Taking a p = 0.51 yields 

<f>„{ u>) = 0.17a fcp (eu o ) 4/3 ur 7/3 . (4.4) 

The result (4.3) has the same form as proposed by Tennekes (1975), while (4.4) 
appears to be new. 
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i i i i > i i li—i 

0 5 10 15 20 25 30 35 40 

(b) kL, 



FIGURE 4. Space-time spectrum of energy, (a) Contours of the energy spectrum 
£(k,w). k is normalized on Lj, and won« 0 /f]- Contours are in factors of \/l0- (b) 
Moments of the frequency distribution of the energy spectrum. Note particularly the 
variance as and the kurtosis Kg. (c) Relation between the space-time spectrum and 
the space spectrum: £(k, u>), (4.2). Results are from a later simulation. 
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FIGURE 5. Auto-correlations in space and time for isotropic turbulence: (a) 
velocity, (b) pressure. 


Figure 5 shows the longitudinal and transverse correlations of the velocity and 
pressure fields. Note that for the velocity, the auto correlation curve Ru(t) has the 
same form as the longitudinal spatial correlation, iZn(ri), and is positive, at least 
for Ti/ < 2.5. These curves are normalized on an integral length scale derived 
from the energy spectrum (Monin & Yaglom 1971, p. 55), 

h = k~ l E(k)dk j J™ E{k)dk. 

For these experiments ij/27r ~ 0.25, i.e., about a quarter of the box size. It is 
convenient to define length and time scales from these one-dimensional correlations 
by a similar procedure as in many experiments. Because of the finite box size, 
the one-dimensional integral scales are estimated from the value of the space or 
time variable where the correlation is 1/e of the variance. So Rn(r = L \{ ) = 
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r/LI 

Figure 6. Comparison of the unnormalized pressure auto-correlation with the 

Hinze /Batchelor theory: computed, theory of Batchelor (1951). Results 

are from a later simulation. 


(l/e)i?n(r = 0). On this basis our computed longitudinal length and time scales 

4V = 1.25IJ, 

r\ E) ~ 1.5Xi/u 0 - 

This value of T\ E) is greater than the value of 1.15li/u 0 obtained by Kyle Squires 
(private communication) in decaying turbulence at Re\ — 35. Both estimates are 
less than the value of 3 inferred by Snyder & Lumley (1971). 

The pressure time and space correlations R pp {t) and R pp {ri) are plotted in fig- 
ure 5b. Note that # pp (ri )/R pp (0) has a similar form as # 22 ( 7 * 1 , 0,0) in that the 
curves become negative when r/Tj >1.5 for the former and r/L\ > 1.3 for the 
latter. A plausible physical explanation can be based on the typical form of the 
large eddies discussed in §2.2. Note that, as for the velocity, the autocorrelation of 
the pressure fluctuations is positive. An important result is that the mean square 
pressure fluctuation is given by 


are 

and 


P 2 /{pul? ^ 1-0. 

This ratio is higher than the estimates made by Batchelor (1951) for low Reynolds 
number wind tunnel turbulence (0.34), and for high Reynolds number turbulence 
by Hinze (1975) (0.5), and George et al. (1984) (0.42). 

A significant agreement with the theory of Batchelor (1951) is shown in figure 6 
where the theoretical value of R pp {r \ ) (using #n( r i ) from a later simulation having 
lower ani sotropy) is compared with computational results. The agreement is quite 
good except near r = 0. 

The final results, figure 7, were the power spectra E\ to E 4 for the four quantities: 
Vp + (u • V)u, Vp, and For high Reynolds number turbulence, where the 
viscous terms are small, the first and third spectra should be equal and also the 
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Figure 7. Power spectra of the terms in the Navier-Stokes equations. 


second and fourth. Furthermore if the eddies are advected past an observer faster 
than they evolve by nonlinear or viscous processes, then following Tennekes (1975) 
a7 and E 3 (k) should be much greater than ^ an< J E 4 (k) for large k. However the 
computations showed that E 3 and E 4 were of comparable magnitude (within 20%). 
But we did find that the ratio E 3 /E 2 increased from about 1 to 8 as kL 3 increased 
from 8 to 36. This is consistent with the advective contribution to (u • V)u being 
much greater than the local nonlinear contribution (which largely determines Vp). 
We also found that E 4 (i.e. was about 2 or 3 times as great as E 2 (Vp), even 

for kL i ~ 20, showing that the contribution of viscous stresses to E 4 and is 
dominant. 

These viscous stresses are probably the explanation for why E 3 is of the same 
order as E 4 . But if the viscous stresses are so relatively large, is this consistent 
with the small eddies being apparently advected by the large eddies? A possible 
explanation is that the small eddies are like slowly decaying vortices advected by 
the large eddies. Within the small eddies the viscous term uV 2 u may be greater 
than the advective acceleration term (U • V)u. But the time for the decay of say an 
extended line vortex with core size k ~ 1 is much longer than (vl l* 2 ) -1 . Therefore we 
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draw the surprising conclusion that, for turbulence in this Reynolds number range, 
the locations of the viscous stresses is associated with relatively large values of -g®- 
suggesting that first there is not necessarily a local instantaneous balance between 
the large straining and dissipation, and second the small eddies decay as they are 
adverted by the large eddies. Perhaps the generation of small scale turbulence 
occurs on shorter time scales than its decay. 

5. Conclusions and Further Work 

1. The central hypothesis to be tested in this project of large scale eddies (or 
Whorls) advecting (or carrying) small scale eddies (or whorls), was supported, if 
not completely confirmed. [Perhaps we should change one word in L.F. Richardson’s 
rhyme: replace ‘Big whorls have little whorls ...’ with ‘... carry ... ’!] 

2. The computations raised numerous questions. To clarify the advection process 
we need to compute the probability distribution of the large eddies, and for the 
pressure correlations the probability distribution of the strain rates. We found that 
the pressure correlation agrees fairly well with the Hinze/Batchelor theory which is 
based on a quasi-normality assumptions for velocity moments. 

3. In these simulations the turbulence was driven by a steady forcing function F. 
The generality of the results is not clear until we explore further the effect of the 
magnitude and spectral distribution of F. Also, do the initial conditions matter? 
However there are encouraging signs of agreement with other kinds of simulations 
and with measurements. These comparisons need to be pursued. 

4. The computations have given rise to some important new physical ideas; the 
connection between the advection of small eddies and the advection of small scale 
pressure fluctuations; the structure of the large scale pressure field; and, finally, 
perhaps some insight, into the temporal pattern of the movement and decay of 
small scale eddies. 

5. There will be several applications of these computations to other areas of tur- 
bulence research, including practical applications. It will be possible to compute 
the effects of turbulence on the generation of internal and surface waves, and to 
develop representations of the space-time structure of turbulence, either as a 3-D 
stochastic field, or by means of conditional eddy techniques. Such representations 
are being used to compute the motions of particles in turbulent flows, where it is 
too expensive to use direct simulations. 
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Table 1 

Numerical and physical parameters 

Number of mesh points: 64 3 

Number of time samples: 288 

Total disk storage used: 604 Mbytes 

Boxsize: 27r 

Force field: F' = J |f| 2 dxj = 0.69 

Kinematic viscosity: u = .02 

Properties of flow field 

Velocity field: = 1.94 

“o = (j«-) 1/2 = 1-14 

Lengthscale: U = t /o°° k ~ lE (k)dk/J™ E{k)dk = 1.25 

Taylor microscale: A = [5 / 0 °° E(k)dk/ J^° k 2 E(k)dkY = 0.84 

Kolmogorov microscale: ^Kol = (i^/e) 1 ' 4 = 0.064 

Reynolds number: Re\ = u 0 \/v = 48 

Dissipation rate: e/(ul/Li) = (15/fJex)(Ti/A) = 0.47 

3 /r 1 3 / 2 

Skewness: | E (fg)*/ (S^) 2 = -0.50 

Normalized force: F'L\ /u\ = 0.66 

‘Integral’ length scales estimated from 1-D correlations: 

Velocity: L\\^/Lj — 1.25, / Lj = 0.76 

Pressure: L ( pp / L\ = 0.75 

‘Integral’ time scales: t\ E ^ /(Li/uq) = 1.5 

T ( 2 E) /{L 1 /«o)= 1-3 
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